A Higher- Accuracy van der Waals Density Functional 
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We propose a second version of the van der Waals density functional (vdW-DF2) of Dion et al. 
[Phys. Rev. Lett. 92, 246401 (2004)], employing a more accurate semilocal exchange functional and 
the use of a large- A'' asymptote gradient correction in determining the vdW kernel. The predicted 
binding energy, equilibrium separation, and potential-energy curve shape are close to those of ac- 
curate quantum chemical calculations on 22 duplexes. We anticipate the enabling of chemically 
accurate calculations in sparse materials of importance for condensed-matter, surface, chemical, 
and biological physics. 



The van der Waals (vdW) attraction is a quantum- 
mechanical phenomenon with charge fluctuations in one 
part of an atomic system that are electrodynamically cor- 
related with charge fluctuations in another. The vdW 
force at one point thus depends on charge events at an- 
other region and is a truly nonlocal correlation effect. 

Methods for accurately calculating the vdW interac- 
tions are critical to understanding sparse matter, includ- 
ing bulk solids {e.g., layered materials, molecular crys- 
tals, and polymers), surface phenomena {e.g., adsorp- 
tion, water overlayers, and gas separation and storage), 
and biostructures {e.g., DNA and protein structure). 

The exact density functional contains the vdW forces. 
Unfortunately, we do not have access to it, but ap- 
proximate versions are abundant. Commonly, the local- 
density approximation (LDA) and generalized gradient 
approximation (GGA) are used with quite some success 
for dense matter, including hard materials and covalently 
bound molecules. They depend on the density in local 
and semilocal ways, respectively, however, and give no 
account of the fully nonlocal vdW interaction. 

First-principles approaches for how vdW can be 
treated in DFT were first proposed for the asymptotic in- 
teraction between fragments jir— These ultimately evolved 
into the van der Waals density functional (vdW-DF) for 
arbitrary geometries.^"*' Despite its success for describing 
dispersion in a breadth of systems better than any other 
nonempirical method^ vdW-DF overestimates equilib- 
rium separations — i^ilr— and underestimates hydrogen- 
bond strength ] ^^'^^ 

In this Letter, we propose a second version of the van 
der Waals density functional (vdW-DF2) employing a 
more accurate semilocal exchange functional PW86i^ii^ 
and the use of a large- asymptote gradient correctionii 
in determining the vdW kernel. By making a full com- 
parison of potential energy curves (PECs) with accurate 
quantum chemistry results for 22 molecular duplexes, we 
show that vdW-DF2 substantially improves (i) equilib- 
rium separations, (ii) hydrogen bond strengths, and (iii) 
vdW attractions at intermediate separations longer than 
the equilibrium ones. The improvement in (iii), found 
via a full PEC comparison, is most critical for important 



"real-life" applications to sparse matter and biological 
matter where it is impossible for basic structural units 
to assume the same separations they would have as bi- 
nary units in vacuo. 

First, we replace revPBE exchange functionaliS with 
PW86,^^'^^ because revPBE is generally too repulsive 
near the equilibrium separation,^ and can bind spuriously 
by exchange alone, although less so than most other lo- 
cal or semilocal functionals. Hence, other exchange func- 
tionals — i2i have been proposed. Recent performance 
studies of various exchange functionals for weakly inter- 
acting atoms2^ and moleculeSfi^ however, show PW86, 
with an enhancement factor proportional to s^/^ at large 
reduced density gradient s, to give the most consistent 
agreement with Hartree-Fock (HF) results, without spu- 
rious exchange binding. It also is a good match^'^ for the 
vdW-DF2 correlation kernel, introduced below, although 
no others were tried. 

The key to the vdW-DF method is the inclusion of 
a long range piece of the correlation energy, i?"'[n], a 
fully nonlocal functional of the density n. This piece is 
evaluated using a "plasmon" pole approximation for the 
inverse dielectric function, which satisfies known conser- 
vation laws, limits, sum rules, and invariancesi^ A sin- 
gle parameter model for the pole position was adopted, 
with the pole residue set by the law of charge-current 
continuity (/-sum rule), and the pole position at large 
wave vector set by the constraint that there be no self- 
Coulomb interaction. The single parameter is determined 
locally from electron-gas energy input via gradient cor- 
rected LDAA 

The nonlocal piece of the correlation energy in both 
vdW-DF and vdW-DF2 is of the form 



= d^r dVn(r)0(r,r')n(r'). 



(1) 



The kernel (f> is given as a function of Rf{r) and Rf{v'), 
where R=|r — r'| and /(r) is a function of 7i(r) and its 
gradient. In fact /(r) is proportional to the exchange- 
correlation energy density exc of a gradient corrected 
LDA at the point r. This arose from the approximately 
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FIG. 1. (Color online) Potential energy curves (PECs) for the best and the worst case of (a) hydrogen-bonded, (b) dispersion- 
dominated, and (c) mixed duplexes. CCSD(T) QC PECs (dash-dotted lines with circles taken from Ref. [24h and the reference 
energies (cross marks taken from Ref. [2^ at the geometry of Ref. [2^ are also shown. The shapes near minima are compared 
in inset figures where PECs are aligned to have the common minimum point. For all the other S22 duplexes, see Supplemental 
Material. 



implemented requirement that the dielectric function im- 
plied by the plasmon pole model should give an exchange- 
correlation energy semilocally consistent with a gradient 
corrected electron gas. We call the semilocal functional 
that fixes /(r) in Eq. ([T]) the internal functional. 

The above is easier to understand for two separate 
molecules, although the arguments apply equally well to 
a pair of high density regions of a sparse material. The 
long range vdW attraction implied by Eq. ([1]) occurs from 
the contribution where r is on one molecule and r' on the 
other. The definition of /(r) and /(r') varies continu- 
ously and independently at each point according to excii^) 
and ea;c(r'). The quantity exc is taken to consist of a gra- 
dient corrected LDA. In the first version of vdW-DF,'* 
the gradient correction was obtained from a gradient ex- 
pansion^ for the slowly varying electron gasi ^^'^^ More 
appropriate is a functional that gives accurate energies 
for molecules, however. When r and r' are each in a sep- 
arate molecule-like region, with exponentially decaying 
tails between them, /(r) and /(r') can both be large and 
give key contributions to a vdW attraction. For this case 
(including perhaps even a molecule near a surface) the 
large- asymptote'^*''^^ and the exchange energy asymp- 
totic series for neutral atoms provide a more accurate ap- 
proximation. In fact, the exchange parameter^ /? of the 
B88 exchange functional, '^^ successfully used in the vast 
majority of DFT calculations on molecules, can be de- 
rived from first principles using the large- A'^ asymptote;^ 
as can the LDA exchange. It seems obvious, then, that 
vdW-DF results should be improved if the second order 
expansion of the exchange in gradients is replaced by the 
second order large- A^ expansion. Interestingly, PW86R 
functional, selected as the overall exchange functional for 
different reasons, also follows the large-A^ behavior for 
small reduced gradient s values down to ^^0.1, where it 
reverts to the form of slowly varying electron gas limit. 

Thus we use 2.222 times larger exchange gradient 



coefficient, a value based on agreement between de- 
rivedii^i^ and empirical^ criteria (a 6% smaller de- 
rived value of Ref. |l7| only gives a marginal improve- 
ment). Assuming that the screened exchange term^^ in- 
creases in the same proportion as gradient exchange it- 
self, finally we get the appropriate gradient coefficient in 
the "Z" notations^ which is multiplied by 2.222. Summa- 
rizing: while Zab = -0.8491 in vdW-DF, Zab = -1.887 
in vdW-DF2, implying changes in the internal functional. 

The performance of our new energy functional is as- 
sessed via comparisons with the accurate S22 reference 
dataset^!^^ and PECs— based on quantum chemistry 
(QC) calculations at the level of CCSD(T) with extrapo- 
lation to the complete basis set limit. These twenty- two 
small molecular duplexes for the non-covalent interac- 
tions typical in biological molecules include hydrogen- 
bonded, dispersion-dominated, and mixed duplexes. Re- 
cent evaluation^ of the performance of the vdW-DF for 
S22 shows it to be quite good, except for H-bonded du- 
plexes, where vdW-DF underestimates the binding en- 
ergy by about 15%. 

Calculations are performed by a plane-wave code and 
an efficient vdW algorithro^i with TrouUier-Martins type 
norm-conserving pseudopotentials. Spot comparison 
with all electron calculations using large basis sets in- 
dicates a calculational accuracy of ~1%, actually better 
than that of most of PAW potentials supplied in vari- 
ous standard codes. Large box sizes were used to control 
spurious electrostatic interactions between replicas. See 
Supplemental Material for further details. 

Figure [T] shows a typical PEC for each kind of inter- 
action. To make a direct comparison to the QC results, 
the PECs are calculated at the same geometry of those 
of the CCSD(T) PEC calculational (shown as dash- 
dotted lines with circles), where each molecule is kept 
in its 822 geometry^^ and moved along the line defined 
by the center-of-mass coordinates of two molecules with- 
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out any rotation. Overall, the vdW-DF2 PECs (solid 
lines) show a remarkable agreement with QC ones for all 
separations and for all three types of interactions. The 
shapes of the PECs near the minima, important for vibra- 
tional frequencies, are greatly improved (see inset figures 
where PECs are aligned to have the common minimum 
point). More importantly the strength of vdW attrac- 
tion at larger distances is weakened in agreement with 
QC, especially for vdW-bonded duplexes (Fig. [TJb)). In 
other words, the original vdW-DF overestimates vdW 
attraction at intermediate separations beyond the equi- 
librium separation, although its minimum energy is ac- 
curate. This has special importance in sparse condensed 
matter. An example will be given below. 

The PECs tend to turn up slightly earlier when ap- 
proaching the repulsion regime at small separations. This 
quite universal feature might be due to PW86 exchange 
functional which is slightly more repulsive than HE at 
short distances. Eor duplexes whose large distance 
asymptote is dominated by dispersion (here methane 
dimer, ethene dimer, and benzene-methane), vdW-DE2 
will have weaker attraction in the asymptotic region and 
smaller Cq coefficients than vdW-DE, which (at least for 
the methane dimer) already gives a Cg coefficient close to 
experiment.® Recently, this prediction was verified in de- 
tail for various duplexes.'^- However, in the region of our 
calculations, such deterioration is not pervasive. In any 
case, neither the vdW-DE curves nor the QC curve have 
reached their asymptotic limit for any of the above three 
cases. The remainder of the 22 duplexes have asymptotic 
forms dominated by electrostatics. 

In Fig. [21 we summarize the calculated binding ener- 
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FIG. 2. (Color online) Comparison for the S22 duplexes of 
(a) binding energy predicted by Jurecka et al.j^S vdW-DF,i 
and the present work (vdW-DF2). (b) binding energy at a 
separation 1 A larger than the equilibrium one and, (c) equi- 
librium separations. The ordinates give the respective devia- 
tions of these quantities from the reference values taken from 
Takatani et al.— for panel (a), Molnar et al.-Si for panel (b), 
and Jurecka et al.— for panel (c). 



gies at equilibrium separation, at equilibrium separation 
plus 1 A, and the equilibrium separations themselves, 
with each quantity given as a deviation from that im- 
plied by the reference calculations^ir— Each subfigure 
clearly shows one of the three major improvements: (i) 
Hydrogen-bond strengths. The mean absolute deviation 
(MAD) of binding energy for hydrogen-bonded duplexes 
is decreased from 91 to 40 meV. (ii) Interaction at inter- 
mediate separations. The MAD of binding energy at a 
separation 1 A larger than the equilibrium one is reduced 
from 35 to 8 meV for the dispersion dominated group and 
also substantially for the mixed group. (Hi) Equilibrium 
separations. MAD is reduced from 0.23 to 0.13 A. Over- 
all, the vdW-DF2 binding energies are within 50 meV of 
the reference except for duplexes 5 and 7. The MAD of 
binding energy is decreased from 41 to 22 meV (13% to 
7.6%). As a final tidbit, we note that in vdW-DF2 the 
MAD for the equilibrium energy of the dispersion dom- 
inated complexes has been reduced to 11 meV, which is 
equal to the MAD of the work of Jurecka et al.r^ which 
until very recently2£ was considered the 'gold standard' 
for accuracy in quantum chemical calculations on this 
group. 

The higher accuracy of vdW-DF2 persists in the ex- 
tended systems we tested. As such applications, we 
calculated: (i) The graphite interlayer binding energy 
and spacing. The binding energy is on a par with 
that of vdW-DF which is already good. Interlayer 
spacing is about 2% shorter than that by vdW-DF, in 
better agreement with experiment, (ii) H2 adsorption 
within two different metal-organic frameworks (MOFs), 
Zn2(BDC)2(TED) and MOF-74. In the former, the bind- 
ing at the highest binding site (57 mcV including zero 
point) is in good agreement with heat of adsorption mea- 
surement (52 meV), whereas vdW-DF overestimates it 
by 60%. This demonstrates the importance of accu- 
rate intermediate-range interaction. For comparison, the 
binding energy of hydrogen to benzene, one of building 
blocks of Zn2(BDC)2(TED), is almost the same for both 
functionals at the equilibrium separation. At larger sep- 
arations, however, vdW-DF substantially overestimates 
the binding. In case of MOF-74, H2 positions are exper- 
imentally known. VdW-DF2 shows its higher accuracy 
in predicting separations between hydrogens and nearby 
atoms (improved by 0.1-0.2 A). The strongest binding 
comes from an unsaturated metal atom, rather than more 
distant structures, and the binding energy, already ac- 
curate in vdW-DF, is little changed in vdW-DF2. We 
also tested vdW-DF2 on the internal structure of a wa- 
ter molecule, as an example of strong chemical bonds in 
molecules. We find that vdW-DF2 is on par with PBE. 
For details of all tests, see Supplemental Material. 

In summary, we have presented an enhanced version 
of vdW-DF, denoted by vdW-DF2, which can be imple- 
mented via simple modifications of existing codes. It re- 
sults in significant improvements in equilibrium spacings 
between noncovalently bound complexes, as well as in 
binding energy, especially when hydrogen bonding plays 
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a role. We make a full comparison of PECs of both 
functionals with accurately known results for a set of 22 
complexes and also apply it to extended solid systems, 
finding favorable results for the new functional. Thus, 
we expect our method to have important applications in 
a wide range of fields, including condensed matter and 
materials physics, chemical physics, and the physics of 



biological materials. 
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Supplemental Material 

Computational details The convergence control parame- 
ters for the plane-wave pseudopotential calculations are 
tuned until the binding energy is converged up to 0.4 
meV (0.01 kcal mol^^). We used 50 Rydberg kinetic- 
energy cutoff for hydrocarbon systems and 60 Rydberg 
for all the others containing nitrogen and oxygen. For 
polar molecules with large dipole moments a 60x60x60 
Bohr'^ cubic unit cell is needed to eliminate spurious elec- 
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trostatic interaction between supercell images. The ac- 
curacy of the pseudopotentials are tested within PBE [J. 
Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 
77, 3865 (1996)] by comparing with all-electron results 
for the water dimer (duplex ^2 in S22) and formic acid 
dimer (duplex #3). The all-electron energies obtained 
were also converged to the level of 0.4 meV (See TableHl). 
For water dimer we obtained binding energies (in meV) 
of 215 (all electron) and 219 (pseudopotential). Similarly, 
for the formic acid dimer we obtained 790 (all electron) 
and 798 (pseudopotential) . The respective deviations are 
thus 2% and 1%. 

Compared to GGA, the vdW-DF2 computation time 
per iteration roughly doubled for the smaller S22 duplex 
calculations, but the increase becomes virtually unmea- 
surable for the larger MOF systems that have more than 
100 atoms per unit cell. The number of iterations re- 
quired for convergence is also typically comparable to 
that of GGA, but in some MOF cases it is almost dou- 
bled. 



TABLE I. Convergence of all-electron PBE binding energy (in 
meV) of the water dimer (duplex #2 in S22) and the formic 
acid dimer (duplex #3). The calculations are performed by 
GaussianOS." Dunning's correlation consistent basis sets with 
diffuse functions are used, with counterpoise-corrections. 



basis set 



water dimer 



formic acid dimer 



aug-cc-pVDZ 
aug-cc-pVTZ 
aug-cc-pVQZ 
aug-cc-pV5Z 
aug-cc-pV6Z 



212 
213 
214 
215 

215 



774 
784 
790 
790 



Gaussian 03, Revision E.Ol, M. J. Frisch, G. W. Trucks, H. B. 
Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, J. A. 
Montgomery, Jr., T. Vreven, K. N. Kudin, J. C. Burant, J. M. 
Millam, S. S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. 
Cossi, G. Scalmani, N. Rega, G. A. Petersson, H. Nakatsuji, M. 
Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. 
Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, M. Klene, 
X. Li, J. E. Knox, H. P. Hratchian, J. B. Cross, V. Bakken, C. 
Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. 
Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, 
P. Y. Ayala, K. Morokuma, G. A. Voth, P. Salvador, J. J. 
Dannenberg, V. G. Zakrzewski, S. Dapprich, A. D. Daniels, M. 
C. Strain, O. Parkas, D. K. Malick, A. D. Rabuck, K. 
Raghavachari, J. B. Foresman, J. V. Ortiz, Q. Cui, A. G. 
Baboul, S. Clifford, J. Cioslowski, B. B. Stefanov, G. Liu, A. 
Liashenko, P. Piskorz, L Komaromi, R,. L. Martin, D. J. Fox, T. 
Keith, M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, M. 
Challacombe, P. M. W. Gill, B. Johnson, W. Chen, M. W. 
Wong, C. Gonzalez, and J. A. Pople, Gaussian, Inc., 
Wallingford CT, 2004. 



TABLE II. Internal structure of a free water molecule. The 
bond length and bond angle calculated by vdW-DF2, vdW- 
DF, and PBE well agree each other within 0.002 A and 0.7 de- 
gree. Norm-conserving pseudopotentials (NCPP), 20x20x20 
Bohr'' cubic unit cell, 130 Rydberg kinetic-energy cutoff, and 
0.002 eV/A force tolerance are used. For comparison all- 
electron (AE) PBE, CCSD(T), and experiment values are 
given as well. AE PBE bond length is 0.007-0.009 A longer 
than our pseudopotential calculation. 



Method 



d(OH) (A) 



6i(H0H) 



vdW-DF2 (NCPP) 
vdW-DF (NCPP) 
PBE (NCPP) 
PBE (AE) 
CCSD(t£ 
Expt."* 



0.960 
0.960 
0.962 

0.969^, 0.971^ 
0.958~ 

0.958 



105.0 

104.7 

104.3 

104.1^ 

104.5" 

104.5 



All-electron calculation with the def2-TZVPP basis set 
(triple-^ quality including high exponent polarization 
functions), which is a larger basis set than the cc-pCVTZ basis 
set, taken from J. Chem. Phys. 126, 124115 (2007). 

All-electron calculation with aug-cc-pVTZ basis set taken from 
J. Phys. Chem. A 108, 2305 (2004). 

D. Feller and K. A. Peterson, J. Chem. Phys. 131, 154306 
(2009). 

S. V. Shirin et al, J. Mol. Spectrosc. 236, 216 (2006). 
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TABLE III. Comparison of binding energy and equilibrium separation. The latter is given as a deviation from the original S22 
geometry" along the center-of-mass line, as shown in the first data column of Table HV)) . 



II 

# 


Duplex 




Binding ener; 


yy (meV) 




Scparat: 


ion (A) 


vdW-DF 


vdW-DF2 


QC^ 




vdW-DF 


vdW-DF2 


1 


Ammonia dimer 


115 


134 


137 


137 


0.20 


0.08 


2 


Water dimer 


185 


218 


218 


218 


0.13 


0.05 


3 


Formic acid dimer 


690 


766 


807 


815 


0.08 


0.08 


4 


Formamide dimer 


587 


655 


692 


699 


0.10 


0.07 


5 


Uracil dimer 


767 


832 


895 


897 


0.07 


0.06 


6 


2-pyridoxine-2-aminopyridine 


639 


687 


725 


737 


0.07 


0.06 


7 


Adenine-thymine 


609 


660 


710 


726 


0.07 


0.05 


8 


Methane dimer 


36 


30 


23 


23 


0.18 


-0.02 


9 


Ethene dimer 


64 


65 


65 


65 


0.36 


0.14 


10 


Benzene-methane 


68 


63 


65 


63 


0.36 


0.20 


11 


Benzene dimer (slip-parallel) 


136 


123 


118 


114 


0.43 


0.30 


12 


Pyrazine dimer 


185 


177 


192 


182 


0.36 


0.24 


13 


Uracil dimer (stacked) 


403 


402 


439 


422 


0.22 


0.11 


14 


Indole-benzene (stacked) 


206 


197 


226 


199 


0.42 


0.29 


15 


Adenine-thymine (stacked) 


461 


466 


530 


506 


0.30 


0.18 


16 


Ethene-ethine 


69 


70 


66 


65 


0.15 


0.03 


17 


Benzene- water 


124 


129 


142 


143 


0.20 


0.06 


18 


Benzene-ammonia 


94 


92 


102 


101 


0.27 


0.13 


19 


Benzene-HCN 


166 


170 


193 


197 


0.24 


0.13 


20 


Benzene dimer (T-shape) 


113 


105 


119 


118 


0.34 


0.22 


21 


Indole-benzene (T-shape) 


214 


206 


248 


243 


0.28 


0.20 


22 


Phenol dimer 


254 


279 


306 


307 


0.19 


0.09 




Mean deviation (MD) 


-36 


-21 


2 




0.23 


0.13 




Mean absolute deviation (MAD) 


41 


22 


7 




0.23 


0.13 




MAD% 


13% 


8% 


2% 









p. Jurecka, J. Sponer, J. Cerny, and P. Hobza, Phys. Chem. Chem. Phys. 8, 1985 (2006). 

T. Takatani, E. G. Hohenstein, M. Malagoli, M. S. Marshall, and C. D. Sherrill, J. Chem. Phys. 132, 144104 (2010). 
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TABLE IV. Left two data columns: The center-of-mass separations and the nearest neighbor (NN) pair separations in the S22 
geometry of Jurecka et al. These center-of-mass distances corresponds to the zeros of the abscissa in the PEC plots in Figs. [TJSl 
The species of NN atom pair are given in the parentheses. Right two data columns: The deviations from the CCSD(T) value 
[J. Chem. Phys. 131, 065102 (2009)] of the binding energy at a distance 1 A larger than that of the S22 equilibrium geometry 
of Jurecka et al. [Phys. Chem. Chem. Phys. 8, 1985 (2006)]. 



# 


Duplex 


Separation 


(A) 




ABE at -1-1 I 


\ (meV) 






Center-of-mass 


NN 


pair 


vdW-DF 


vdW-DF2 


1 


Ammonia dimer 


3.21 


2.50 


(NH) 


10 


3 


2 


Water dimer 


2.91 


1.95 


(OH) 


8 


7 


3 


Formic acid dimer 


2.99 


1.67 


(OH) 


6 


16 


4 


Formamide dimer 


3.23 


1.84 


(OH) 


5 


7 


5 


Uracil dimer 


6.07 


1.77 


(OH) 


-4 


-5 


6 


2-py ridoxine • 2- aminopy ridine 


5.14 


1.86 


(NH) 


6 


3 


7 


Adenine -thymine 


5.97 


1.82 


(NH) 


N/A^ 


N/A^ 


8 


Methane dimer 


3.72 


3.51 


(CH) 


13 


2 


9 


Ethene dimer 


3.72 


2.56 


(HH) 


21 


4 


10 


Benzene-methane 


3.72 


2.79 


(CH) 


19 


4 


11 


Benzene dimer (slip-parallel) 


3.76 


3.37 


(CC) 


40 


12 


12 


Pyrazine dimer 


3.48 


3.27 


(NH) 


42 


10 


13 


Uracil dimer (stacked) 


3.17 


2.71 


(OH) 


63 


12 


14 


Indole-benzene (stacked) 


3.50 


3.20 


(CH) 


50 


10 


15 


Adenine-thymine (stacked) 


3.17 


2.68 


(HH) 


N/A^ 


N/A^ 


16 


Ethene-ethine 


4.42 


2.83 


(CH) 


12 


3 


17 


Benzene-water 


3.38 


2.60 


(CH) 


13 





18 


Benzene-ammonia 


3.56 


2.77 


(CH) 


16 


1 


19 


Benzene-HCN 


3.95 


2.67 


(CH) 


12 


-5 


20 


Benzene dimer (T-shape) 


4.91 


2.80 


(CH) 


20 


-3 


21 


Indole-benzene (T-shape) 


4.88 


2.59 


(CH) 


10 


-19 


22 


Phenol dimer 


4.92 


1.94 


(OH) 


N/A^ 


N/A^ 




Mean deviation (MD) 








19 


3 




Mean absolute deviation (MAD) 








19 


7 




MAD% 








16% 


4% 



^ The CCSD(T) PEC is not available. 



TABLE V. Graphite interlayer binding energy and interlayer spacing (i.e., half of the c lattice constant). An AB stacked 4-layer 
unit cell is used for vdW-DF and vdW-DF2 calculations. The interlayer binding energy is calculated by subtracting the energy 
of each single layer in the same unit cell from that of 4-layer infinite bulk. 





vdW-DF 


vdW-DF2 


QMC^ 


Experiments^ 


Binding energy (meV/atom) 


50 


49 


56±5 


52±5 


Interlayer spacing (A) 


3.60 


3.53 


3.43 


3.36 



L. Spanu, S. Sorella, and G. Galli, Phys. Rev. Lett. 103, 196401 (2009). Zero-point energy and phonon contributions at 300K are 
included. 

R. Zacharia, H. Ulbricht, and T. Hertel, Phys. Rev. B 69, 155406 (2004). 
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FIG. 1. Potential energy curves (PECs) of hydrogen-bonded duplexes. The S22 data points (cross marks) are taken from 
T. Takatani, E. G. Hohenstein, M. MalagoU, M. S. Marshall, and C. D. Sherrill, J. Chem. Phys. 132, 144104 (2010). The 
CCSD(T) PECs data (dashed lines with open circles) in this and subsequent figures are taken from L. F. Molnar, X. He, 
B. Wang, and J. Kenneth M. Merz, J. Chem. Phys. 131, 065102 (2009). For the hydrogen-bonded adenine-thymine duplex, 
CCSD(T) PEC data is not available. The shapes near minima are compared in inset figures where PECs are aligned to have 
the common minimum point. 
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FIG. 2. Potential energy curves (PECs) of dispersion-dominated duplexes. For the stacked adenine-thymine duplex, CCSD(T) 
PEC is not available. The S22 data points (cross marks) are taken from T. Takatani, E. G. Hohenstein, M. Malagoli, M. S. Mar- 
shall, and C. D. Sherrill, J. Chem. Phys. 132, 144104 (2010). The original S22 values (plus marks taken from P. Jurecka, 
J. Sponer, J. Cerny, and P. Hobza, Phys. Chem. Chem. Phys. 8, 1985 (2006)) are also shown if the difference between those 
two values are larger than 3%. 
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FIG. 3. Potential energy curves (PECs) of mixed interaction duplexes. 
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FIG. 4. Comparison of vdW-DF with revPBE exchange functional (vdW-DF/revPBE), vdW-DF with PW86R (vdW- 
DF/PW86R), and vdW-DF2 with PW86R (vdW-DF2). Potential energy curves (PECs) of water dimer (hydrogen-bonded du- 
plexes), ethene dimer (dispersion), ethene-ethine duplex (mixed) are shown. Note that the difference between vdW-DF/revPBE 
and vdW-DF/PW86R vanishes beyond 1 A. Replacing revPBE exchange with PW86 weakens the repulsion near equilibrium 
separations and gives more accurate separations but also large overbinding. Furthermore, the comparison of PECs with 
CCSD(T) shows that the original vdW-DF consistently overestimates the vdW attraction. The use of large-N asymptote in 
determining the internal functional for vdW kernel effectively weakens the vdW attraction and vdW-DF2 gives an excellent 
agreement with CCSD(T) for all separations. 
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Center-to-center Distance (A) 



FIG. 5. Hydrogen-benzene potential energy curve (PEC). Sampled along a normal vector of the benzene plane passing through 
the benzene center. H2 is parallel to the vector. The CCSD(T) result is taken from Hiibner et al. J. Phys. Chem. A 108, 3019 
(2004) and QMC from Beaudet et al. J. Chem. Phys. 129, 164711 (2008). 




FIG. 6. Hydrogen adsorption by the metal-organic framework (MOF) structure Zn2(BDC)2(TED) where BDC and TED denote 
the linking structures so labeled above. Hydrogen atoms attached to carbon atoms are removed for simplicity. The left panel 
shows the H2 in a "corner" site along the body diagonal, and the right panel shows the H2 head-on in the "face center" site. 
The binding of the H2 is dominated in both cases by its interactions with the benzene and TED complexes which are ~5 A 
away. 
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TABLE VI. Hydrogen binding energies in Zn2(BDC)2(TED) 
(in meV). Zero-point energies are included in the calculated 
values. These values demonstrate the importance for the func- 
tional to give accurate interaction energies at intermediate 
distances that are longer than the equilibrium distances in 
isolated duplexes. For example, the hydrogen-benzene inter- 
action (Fig.[SJ at a center-to-center distance of 5 A is ~100% 
overestimated by vdW-DF, although the binding energy at 
the equilibrium separation is quite accurate. This effect leads 
to the 60% overestimate of the binding of the most strongly 
bound H2 in the Zn2(BDC)2(TED) crystal by vdW-DF. In 
vdW-DF2 the error reduced to a reasonable size, and a switch- 
ing of which site is the most strongly bound also occurs. 

Ho hiudin? site H2 binding energy (meV) 

H2bmdmgsite ^^^.ppa ^dW-DF2 Expt.^ 

Corner site 74 65 

Face-center site 82 56 

^ L. Kong, V. R. Cooper, N. Nijem, K. Li, J. Li, Y. J. Chabal, 
and D. C. Langreth, Phys. Rev. B 79, 081407{R) (2009). 

^ J. Y. Lee, D. H. Olson, L. Pan, T. J. Emge, and J. Li, Adv. 
Funct. Mater. 17, 1255 (2007). 




FIG. 7. Hydrogen adsorption sites (black dots) in MOF-74 [N. L. Rosi, J. Kim, M. Eddaoudi, B. Chen, M. O'Keeffe, and O. 
M. Yaghi, J. Am. Chem. Soc. 127, 1504 (2005)]. (Left) Top view into the pore. (Right) Side view of the pore waU. The 
viewing direction and area are depicted by a short arrow and a rectangle, respectively, in the left figure. 
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TABLE VII. Equilibrium separations (in A) of hydrogens to the nearest atom, to the nearest hydrogen, and to nearby benzene 
(Bz) rings in MOF-74. Calculated binding energies (including zero-point energies) of H2 at Zn site are given at the last line. 
In this MOF, the binding is dominated by the proximity of the H2 to the unsaturated Zn atom, and the accurate binding 
energy predicted by vdW-DF is essentially unchanged in vdW-DF2. Still, the weakening of the intermediate range attraction 
is presumably compensated by the closer proximity of the H2 to the Zn atom in vdW-DF2. 



From 


To^ 


vdW-DF 


Distance (A) 






vdW-DF2 


Expt._^ 




Zn 


3.0 


2.8 


2.6 


Zn H2 site 


Bzl^ 


4.2 


4.1 


3.9 




Bz33 


4.4 


4.3 


4.2 




H2 Tite 


3.3 


3.1 


2.9 




01 


3.7 


3.4 


3.3 


H2 site 


02 


4.3 


4.0 


3.5 


03 


3.4 


3.3 


3.1 




Bz2^ 


5.1 


4.8 


5.0 




BzS^ 


4.7 


4.7 


4.6 




Bz4~ 


5.4 


5.3 


4.7 


H2 binding 


; energy at Zn site (eV) 


0.10 


0.10 


0.09 



L. Kong, G. Roman-Perez, J. M. Soler, and D. C. Langreth, Phys. Rev. Lett. 103, 096103 (2009). 
^ Y. Liu, H. Kabbour, C. M. Brown, D. A. Neumann, and C. C. Ahn, Langmuir 24, 4772 (2008). 
See the right panel of Fig. [7] for a depiction of the corresponding atomic locations. 
To the center of the molecule. 




FIG. 8. D2 adsorption sites in the fully loaded MOF-74: Via neutron diffraction" (left panel); Via vdW-DF2 (right panel). 
"Y. Liu, H. Kabbour, C. M. Brown, D. A. Neumann, and C. C. Ahn, Langmuir 24, 4772 (2008). 



